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Abstract 

The macroscopic quantum interference effects in ballistic Joseph- 
son microstructures are investigated. The studied system consists 
of four bulk superconductors (terminals) which are weakly coupled 
through the mesoscopic rectangular normal metal (two dimensional 
electron gas) . We show that nonlocal electrodynamics of ballistic sys- 
tems leads to specific current-phase relations for the mesoscopic multi- 
terminal Josephson junction. The nonlocal coupling of supercurrents 
produces the "dragging" effect for phases of the superconducting order 
parameter in different terminals. The composite Josephson junction, 
based on this effect, exhibits the two -level system behaviour controlled 
by the external magnetic flux. The effect of magnetic flux transfer in 
a system of nonlocally coupled superconducting rings is studied. 

Keywords: MultiterminalJosephson junction; Nonlocal coupling; Drag effect; 
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1 Introduction 



Recently there has been considerable interest in ballistic mesoscopic S-2DEG- 
S Josephson junctions (2DEG - two dimensional electron gas) . Heida et.al. 
H] investigating this new class of fully phase coherent Josephson junctions 
of comparable width and length, have demonstrated the nonlocal phase de- 
pendence of mesoscopic supercurrents. Nonlocality leads to the anomalous 
critical current dependence on the magnetic flux in such mesoscopic junctions 
[1-3]. It was noted in ^ that the four-terminal generalization of S-2DEG- 
S structures (namely, considering of the case when two-dimensional normal 
layer is connected with four terminals instead of two ones as in Refs.[l-3]) 
will allow the studying of the specific multichannel macroscopic quantum 
interference effects in mesoscopic multiterminals. For review of experimen- 
tal and theoretical investigation of the conventional multiterminal Josephson 
junctions see P] and references there. The microscopic theory of coherent 
current states in mesoscopic ballistic 4-terminal Josephson junction was of- 
fered in 0]. In present paper we study the macroscopic manifestations of 
the nonlocal coupling of supercurrents inside the mesoscopic multiterminal 
weak link. We show, that 4-terminal Josephson junction with two terminals 
incorporated in superconducting rings, presents the macroscopic bistable sys- 
tem in which the energy levels and the potential barrier height are regulated 
by the magnetic flux $ applied through the ring. Such composite S-N($)- 
S Josephson junction between conventional superconductors behaves similar 
to so called "vr-contact" (see the recent paper [0 on tunable qubit, based 
on "yr-contact" between d-wave superconductors). In addition, we introduce 
the effect of magnetic flux transfer in system of two superconducting rings 
weakly coupled through the mesoscopic 4-terminal junction. 



2 System description. Basic equations 

The mesoscopic four-terminal Josephson junction we study is sketched in Fig. 
la. It consists of four bulk superconductors (S) which are weakly coupled 
through a ballistic rectangular normal (N) layer (two dimensional electron 
gas ) of length L and width W . Contact sizes L, W are suppose to be 
much larger than Fermi wave length A^. In the stationary state, the ballistic 
motion of the electrons through the normal junction region (— L/2 < x < L/2 
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-W/2 < y < W/2) is described by the quasi-classical Eilenberger equations 



d ^ 

y^Fg^G+[iufs + A,G] = 0, 

The electrons go through quasi-classical trajectories which are straight 
lines along the electron velocity on the Fermi surface, v^. The trajectory 
which is coming from the ith bank and going to the jth bank {i j tra- 
jectory) is shown in Fig. lb. In the point p{x, y) , the interference between 
the trajectories connecting different superconducting banks {e.g. i ^ j and 
i' —>■ j', in Fig. lb) takes place. The state of zth S-terminal is defined by the 
complex pair potential A = Aoe*"^' {i = 1, ...,4) ,i.e. by the phase (pi in the 
ith bulk superconductor. 

The solution of Eilenberger equations (1) inside N-layer can be easily 
obtained by integrating over the electronic trajectories in the same way as for 
a ballistic point contact 0. In Ref. 0, the anomalous fuj{'^F,p) and normal 
QujiyPiP) Green's functions were calculated for arbitrary temperatures < 
T < Tc and arbitrary sizes L,W (in the scale ~ ^vp/^o) • For the 
current density j(p) and self-consistent off-diagonal potential A, we have the 
expressions: 

j(p) = -47r^eiV(0)T^(v^(7^), (2) 

(jj>0 

A = A27rT Y: {U) ■ (3) 

Here uj is Matsubara frequency, A^(0) = m/(27r), (...) is the averaging over 
directions of 2D vector vj?, A is the constant of electron-phonon coupling. 
Inside the normal layer A = and hence A = 0. However, the anomalous av- 
erage, {ipip) has the finite value and determines the induced order parameter, 
*(p), as 

vl/(p-)=27rT5:(/^). (4) 
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3 Current-phase relations 



We consider the case of small junction {L,W « ^o) when the effects of 
nonlocality are most pronounced. Qualitatively, it describes the situation of 
L,W ^0- For general consideration see [^. The total current flowing into 
the iih terminal, /j, depends on the phases ipj {j = 1, ...,4) in all the banks: 

vrAo^ . ^i-ipj Aq cos(^^^) 
I^ = 2^ 1^3 sm ( tanh [ ] . (5) 

In the case of two-terminal junction we have from (5) the Josephson 
current of a clean point contact]^. The geometrical factors, 7jj, present the 
coupling coefficients between the iih and jth terminals and are functions of 
the length to width ratio k = L/W as 

e'^ppd, 1 e'^ppd . k 

7i3 = — ^(1-^^^)' 724 = — ^(1 



e^ppd 1 + k 

712 = 714 = 723 = 734 = -^^(-^^== - 1), %j = -fji. (6) 

where d = a/L^ + and the numeration of the banks is shown in Fig. lb. 
In eq. (5), the positive sign of li corresponds to the direction of the current 
from the normal layer to the ith bank. Note that Y^i=i li = 0. 
In the following we consider, for simplicity, the symmetric junction (square 
N-layer), k = 1, and the case of temperatures T near the critical temperature 
Tc. We emphasize, that the studied here macroscopic interference phenomena 
take place at low temperatures, e. g. at T = 0, as well as near Tg. The 
difference is in the values of the critical currents and in the replacing of the 
sinusoidal current-phase dependence sinip by sin((^/2) sign[cos {ip/2)]. 
Near Tc, and for = 1, eq. (5) takes the form 

h = 71" ~ ""'^^ ^ ^ 

with 

7l3 = 724 = V2, 7l2 = 7l4 = 723 = 734 = 1- (S) 

The corresponding Josephson coupling energy of the 4-terminal junction, Ej, 
(such that li = {2e/h)dEj/dipi) has the form 
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2i 4eTc ^2 



j<fe 



^7jfc[l -cos(¥?j -y^fc)]. (9) 



In the case of conventional 4-terminal Josephson junction [^], which presents 
the system of crossed dirty superconducting microbridges between bulk su- 
perconductors, the current-phase relations have the form similar to (7). The 
essential difference consists in the structure of the coupling coefficients •jij. 
In the mesoscopic ballistic junction these form factors are not factorized 
(namely, 712734 7^ 713724)5 in spite of the conventional case, where •jij ~ 
{1/ Ri){l/ Rj); Ri are the normal resistances of the microbridges. 

The current-phase relations (7) determine the behaviour of the system in 
the presence of the transport currents and (or) the diamagnetic currents in- 
duced by the magnetic fluxes through closed superconducting rings. We will 
distinguish two types of the circuit implication of the mesoscopic 4-terminal 
junction (Fig. 2). The first one, "crossed" or "transverse" implication is 
shown in Fig. 2a. It closely corresponds to the case of conventional 4- 
terminal junction. Thus, the theory developed in |^ is applicable to this 
case too. In the "parallel" (Fig. 2b) case, the new specific features in the be- 
haviour of the 4-terminal Josephson junction appear. They are results of the 
nonlocal coupling of currents inside the N-layer, which produces the peculiar 
effect of " dragging" of the phase difference between one pair of terminals by 
the phase difference between another pair of terminals. We will concentrate 
here on the "dragging" case (Fig. 2b). 

By introducing the phase differences 
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and taking into account the relations 



I = l2 = -h 



J = h = -L 



4, 



(10) 



we obtain from (7) 



J = sin 61+ [(^2 + l)sin-cos- + (v^- 1) cos - sin -] cosx, (H) 

2 2 2 2 
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J — sin if + [(v^ + 1) sin — cos - + (v^ — 1) cos — sin -] cos x- (12) 

All the currents are normalized by 

_ 7rAo(r)^eVL(v/2-l) 

Prom the current conservation, it follows that the phase x iii ©qs. (11) and 
(12) can takes only two value of or tt. The dimensionless Josephson coupling 
energy Ej (in units hji^e)!^ in terms of Q and x takes the form 



Ej — — cos^ — cos</7 — [(•\/2 + l)cos — cos- — (-\/2 — l)sin — sin-]cosx. (13) 

Z Z Z Z 

The minimization of Ej with respect to x actually gives that phase x takes 
the value or tt, depending on the equilibrium values of Q and y?: 

cos X — sign\{\Pl + 1) cos ^ cos - — (a/2 — 1) sin ^ sin -]. (14) 

The coupled system of equations (11), (12) with condition (14) for x de- 
termines the current-phase relations l{0](p) and J{ip;9). We note that in 
the "crossed" implication case (as well as in the conventional 4-terminal 
Josephson junction) the coupling terms in eqs. (11)- (14), which contain the 
coefficient (a/2 — 1), are absent. The appearance of this additional couphng 
leads to new effects. 



4 Bistable states in composite S-N($)-S junc- 
tion 

Let us consider the case when the terminals 3 and 4 are incorporated in a 
superconducting ring. When the external flux is applied, the phase ip is con- 
trolled by the magnetic flux which is embraced by the ring: (p — 27r$/$o 
($0 = h/2e). The configuration shown in Fig. 3 presents the Josephson weak 
link (the region enclosed by the dashed rectangular) between two supercon- 
ductors. This weak link has the current-phase relation I (9), which depends 
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on the internal degrees of freedom (trapped flux $) of such composite S- 
N($)-S junction. The nonlocal coupling in the N-layer leads to the nontriv- 
ial behaviour of the I{9) as function of It follows from (11) and (12) 
that for $ 7^ 0, due to the term proportional to (\/2 — 1), zero value of 
the current / corresponds to nonvanishing value of the phase 9. In Fig. 4 
the function I{9) is shown for two values of flux, $ = (dashed curve) and 
$ = $o/2 (solid curve). The phase dependence I (9) is splitted near 9 = 
at $ 7^ 0. If $ = $o/2 and 7 = the phase difference is not zero but 
equals to ^ = ±9o = ± arcsin ( ^~| ). Both values ±9o correspond to the 
equal minima of potential Ej{9\^ = $o/2) as shown in Fig.5. With going 
of $ apart from ^o/2 the potential profile is changed as shown in Fig.5, and 
one of the minima disappears. Thus, the S-N(<l>)-S junction presents the two 
level system with energy levels which can be regulated by the magnetic flux 
At $ = $o/2 we have the bistable system with equal energy values of the 
levels. The picture resembles the so called "tt- contacts", i.e. the conventional 
Josephson junction with magnetic impurities in barrier or contact between 
d-wave superconductors. 



5 Magnetic flux transfer 

The appearance of nonzero phase difference 9 at zero value of current I 
in S-N($)-S system leads to the effect of magnetic flux transfer, which is 
a macroscopic manifestation of the nonlocal coupling in the mesoscopic 4- 
terminal Josephson junction. 

Consider the case when both pairs of terminals 1-2 and 3-4 are short 
circuited by superconducting rings ( S-N($)-S junction incorporated in the 
ring). The currents I and J in eqs. (11), (12) are now related to the corre- 
sponding magnetic fluxes through the rings by 

I-—r , J-—r , (15) 

where L/, Lj are self-inductances of the rings, and are the external 

magnetic fluxes applied to the rings and and $ are the resulting embraced 
fluxes. The phases 9 and ip are related to the magnetic fluxes by ^ = 27r6/$o, 
cp = 27r$/$o- 
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From eqs. (11), (12) and (15) we obtain the coupled equations for 9 and 



e = e'^-/:7{sine+[(V2+l)sin® cos^+(V2-l)cos® sin^]cosx}, (16) 

r- $ e ^ $ e 

$ = $'^_/:^{sin$+[(V2+l)sin-cos- + (V2-l)cos-sin-]cosx}, (17) 

where we have introduced dimensionless self-inductances = {2el?i)lQLi j 
and magnetic fluxes are measured in units h/2e. It follows from these equa- 
tions that, when the flux is applied to the J- ring, even at 0*^ = the 
magnetic flux in /-ring 6 is not zero. Suppose for simplicity that Cj « 1, 
so that $ ~ For the transferred flux (at B*^ = 0) we obtain the 
equation: 

r- 

e = -Ci{sm. e + [(V2 + 1) sin — cos —+ 

, ^ , e ^ , e , ^ , 0,, 

(v 2 — 1) cos — sm — J sign[[-\/2 + 1) cos — cos (v 2 — 1) sin — sm — J|. 

2 2 2 2 2 2^^^ 

In the limit £/ — 0, we obtain from (17) for transferred flux 

e = -/:/(V2 - l)siny sic/n(cosy). (19) 

The numerical solutions of eq. (18) for different values of £/ are presented 
in Fig. 6. The transferred flux B is a periodical function of applied flux <l>'^, 
G(<l'^ + 2n7r) = B($'^) where n is an integer number. For finite values of £/, 
the small hysteresis region exists in the dependence 0($'^) near $^ = $o/2 
(see Fig. 6). 



6 Summary 

The multiterminal Josephson junction presents the system in which weak cou- 
pling takes place between several massive superconducting banks (terminals). 
The general phenomenological theory of Josephson multiterminals |]^ shows 
the variety of macroscopic quantum interference phenomena in such systems. 
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In present paper, we have introduced and studied a new kind of Josephson 
multiterminals, the mesoscopic balhstic multiterminal Josephson junction. 
In the studied system the nonlocal coupling is established through a meso- 
scopic two dimensional normal metal layer. Within the microscopic theory, 
the current-phase relations for mesoscopic multiterminal junction have been 
calculated. They take into account the nonlocal coupling of supercurrents 
inside the ballistic weak link. This current-phase relations describe specific 
(as compared to the conventional multiterminals 0) behaviour of the system 
in the presence of the transport currents and (or) the diamagnetic currents 
induced by the magnetic fluxes through the closed superconducting rings. 
We showed that nonlocality produces the dragging of the phase difference 
between one pair of terminals by the phase difference between another pair 
of terminals. 

We have studied the macroscopic manifestations of this dragging in two 
situations. The first one is, what we call S-N($)-S Josephson junction (Fig. 
3). It presents the S-N-S contact with a N-layer which properties are reg- 
ulated by the external magnetic flux applied through the ring as shown in 
Fig. 3. The internal degree of freedom (magnetic flux $ embraced in the 
ring) makes this macroscopic system similar to the vr-contacts, i. e. to 
the contacts having non-zero ground state phase difference at zero value of 
the current. We obtained that our " macroscopic vr-contact" presents the 
bistable system in which the energy levels and potential barrier are regu- 
lated by the magnetic flux $. At $ equals to 1/2 of magnetic flux quantum 
$0; we have doubly degenerate ground state, just as in the case of vr-contact 
between d-wave superconductors (see 0). 

Secondly, we consider a system of two superconducting rings weakly cou- 
pled through the mesoscopic 4-terminal junction. Such system was studied 
in 1^ for the case of conventional 4-terminal Josephson junction. Now (in 
addition to the results of Ref. 0), in nonlocal mesoscopic case the effect of 
magnetic flux transfer is obtained. It consists of appearance of magnetic flux 
in one of the rings when the magnetic flux is applied through the another 
ring. The amplitude of effect is calculated for different values of the ring's 
selflnductance. 

In conclusion, it is interesting to realize the 4-terminal SQUID device 



Toll based on mesoscopic 4-terminal junction, which will expose the studied 



here effects of nonlocality. 
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Figure captions 

Fig. 1. (a) The mesoscopic ballistic 4-terminal Josephson junction. Four 
bulk superconductors (S) are weakly coupled through a rectangular normal 
(N) layer of length L and width W . b) The electron quasi-classical trajec- 
tories inside the normal region (— L/2 < a; < L/2, —W /2 < y < W/2) are 
straight lines along the electron velocity on the Fermi surface \p. The tra- 
jectory i ^ j connects the ith and jth sides. 

Fig. 2. Two types of implication of mesoscopic 4-terminal Josephson 
junction in the external circuits, (a) Crossed (b) Parallel. 
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Fig. 3. The S-N($)-S Joscphson junction. Two superconductors (S) are 
weakly coupled through the composite weak link (dashed rectangular) which 
consists of the N-layer and a coupled superconducting ring with the embraced 
magnetic flux 

Fig. 4. The current I{9) through the S-N($)-S Joscphson junction is 
shown for two values of $: $ = (dashed curve) and $ = $o/2 (solid 
curve). For $ ^ $o/2, I{0) is splitted near the 6* = 0. 

Fig. 5. The Joscphson coupling energy for a S-N($)-S junction for 
$ = $o/2 (solid curve) has two equals minima. For $ > $o/2 or $ < $o/2 
(dashed curves) the potential profile is changed and one of the minima dis- 
appears. 

Fig. 6. The dependence of the transferred flux to the /-ring, G, on the 

flux in the J-ring, $, for different values of selfinductance Cf. " " - Ci = 1, 

"o o o" - Ci = 0.5, - Ci = 0.1. The dashed rectangular shows the 

smal hysteress region near $ = tt. 
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